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New limb-darkening coefficients and synthetic photometry 
for model-atmosphere grids at Galactic, LMC, and SMC 
abundances. 
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ABSTRACT 

New grids of Atlas9 models have been calculated using revised convection parame- 
ters and updated opacity-distribution functions, for chemical compositions intended 
to be representative of solar, [M/H] = +0.3, +0.5, Large Magellanic Cloud (LMC), 
and Small Magellanic Cloud (SMC) abundances. The grids cover T e g = 3.5-50kK, 
from log g = 5.0 to the effective Eddington limit. Limb-darkening coefficients and 
synthetic photometry are presented in the UBVRIJHKLM , uvby, ugriz, WFCAM, 
Hipparcos/Tycho, and Kepler passbands for these models, and for Castelli's compara- 
ble 'new-ODF' grids. Flux distributions are given for the new models. The sensitivity 
of limb-darkening coefficients to the adopted physics is illustrated. 



1 INTRODUCTION 



The significance of limb darkening as a probe of atmospheric 
structure was recogni zed from the earliest modelling of stel- 
lar at mospheres (e.g. JSchwarz schild 1906; Jeans 1 9171 ; iMilnel 
1 192lh . and its importance in the quantitativ e photometric 
analysis of eclipsing bin aries similarly noted (| Russelll 1 1 9 1 2| ; 
iRussell fc Shaplevlll912l ). Reliable empirical determinations 
of limb darkening are rarely possible in eclipsing binaries (in 
part because of degeneracies with other parameters), and 
so values derived from model-atmosphere calculations re- 
tain their importance through to the present day, where 
new applications include modelling of exoplanetary tran- 
sits, microl ensing events, and spatially re s olved stellar sur- 
faces (e.g.. ISouthworthl |200S|; IWittl 1 19951 : iHestrofferl fl997l ; 
lAufdenberg. Ludwig fc Kervellall2005l ). 

Wholesale calculation of limb- darkening coefficients be- 
came feasible only after production of the first large grids of 
model atmospheres (cf. ICrygar. Cooper fc Jurkevichlll972l . 
who review earlier studies). In particul ar, the exten 



sive model-atmosphere grids ge nerated by Kurucj d 19791 . 



1993T) form the basis of wo r k by Wade fc Rucinski 



Dlaz-Cordoves fc Gimened (| 19921) , Ivan Hammd (1993 
Claret] (|2000h , and othersT 



cheap, to generate self-consistent Atlas grids covering large 
ranges in parameters of interest. The aspect of internal con- 
sistency is of particular importance when modelling the 
spectra of systems whose temperatures and gravities may 
show considerable variations over their surfaces (such as 
the components of close binary systems, rapidly rotating 
stars, and non-radial pulsators), and when attempting uni- 
form analyses of samples encompassing a range of stellar 
parameters (for example, in synoptic studies of light-curves 
of exoplanetary transits). 

Almost all extensive calculations of grids of Atlas- 
based limb-darkening coefficients published to date have 
used the formulations of opacity-distribution functions 
(ODFs) and atmos pheric co n vectio n inherent to the mod- 
els published by iKurucj (|l993h Fl However, improved 
treatements in both areas have emerged subsequently, 
leading in some cases to quite si gnificant changes in 
the models (|Castelli fc Kuruczl [200? ) . The principal pur- 
pose of the present paper is to present calculations 
of limb-darkening coefficients made using new ODFs 
and improved parameterizations of convection. Because 
of growing interest in mod elling eclipsing binaries in 
the Magellanic Clouds (e.g.. | Harries. Hilditch fc Howarthl 



Although several programs are now available in the 
public domain for generating line-blanketed model atmo- 
spheres which relax the approximation of Local Thermo- 
dynamic Equilibrium (LTE), Kurucz's Atlas codes remain 
the benchmark for LTE calculations, and have two major 
attractions for the present purposes. First, while the mod- 
els may not accurately reproduce individual spectral lines 
for which non-LTE effects are significant, they nonetheless 
appear to be quite successful in generating accurate atmo- 
spheric structures and broad-band fluxes. 

Secondly, it is straightforward, and computationally 



20031; iHilditch. Howarth fc Harries! 120051 ; iBonanosI bOOoT 
North et all |2010h . new models with tailored MC abun- 
dances are included. 



2 MODELS 

Atlas has undergone more or less continuous devel- 
opment over its ~40-year lifetime, and currently ex- 

1 An exception is work reported bv lBarban et al" 
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ists in two distinct forms, Atlas9 (which uses opacity- 
distribution functions) and Atlas12 (which uses opacity 
sampling). The models discussed here were all computed 
using the Trieste port of Atlas9 t o GNU-linux systems 
(|Sbordone, Bonifacio fc Castellill2007h . 

2.1 ODFs and Abundances 

Opacity-distributi on functions are required as a basic input 
to Atlas9 models. ICastelli fc Kuruczl (120041 ) describe ODFs 
which incorporate a v ariety of improvements over those used 
in 

theH 

urucj (|l993l ) grids, of which perhaps the most im- 
portant is the inclusion of improved and additional molecu- 
lar opacities. The model s described he r e use new OD Fs, gen- 
erated using dfsynthe dKuru"c3l2005l ; ICastellil 120051 ') , based 
on the lCastelli fc Kurucj (|2004T l line lists. 

The baseline models use solar abundances reported by 
lAsplund. Grevesse fc Sauvall (|2005l ). Limb darkening plays 
an important role in the analysis of light-curves of stars with 
transiting exoplanets, and the sample of such stars currently 
known shows a propensity for enhanc e d metailicity (e.g . , 
lGonzalez|[l997l ; iFischer fc Valentill2005l ; ISousa et al]|200sh . 
With this in mind, additional models (and ODFs) were cal- 
culated with metallicities enhanced by +0.3 and +0.5 dex 
over solar values. For test purposes, ODFs (but not model 
grids) were computed at several other metallicities, including 
[M/H] = —0.5 for the Vega models discussed in Section^ 

The routine study of eclipsing binaries in external galax- 
ies, and especially in the Magellanic Clouds (MCs), has be- 
come feasible in recent years. While no one set of MC abun- 
dances is likely to be generally satisfactory, most interest in 
eclipsing binaries in the Clouds has concentrated on early- 
type stars (in part because of their intrinsic brightness). In 
order to have an appropriate (Pop. I) reference set of rep- 
resentative Large and Small Magellanic Cloud abundances 
that are at least well defined, and apparently reasonable, 
results have been compiled from a number of sources. These 
abundances, listed in Appendix [Al (Table lATj) , have been 
used to compute new LMC and SMC ODFs. 

All the ODFs newly calculated here are available on- 
line. 



2.2 Convection 

All model-atmosphere codes incorporate approximations 
and parameterizations of physics that may be difficult or ex- 
pensive tomodel, or is simply poorly understood (cf., e.g., 
lKuruczlll996l ). Convection is a prime example, and is char- 
acterized in Atlas9 by a mixing-length approximation with 

optional overshooting. 

The widely distributed iKuruc 3 (|l993h grids were 
calculated with a ratio of mixing length to pressure 
scale height of l/H = 1.25 a n d ove rshooting included, 
but ICastelli. Gratton fc Kuruczl (1 19971) argue that more- 
consistent results are obtained with Atlas9 when over- 
shooting is switched off. Furthermore, ISmallevI (|2005l . and 
personal communication) advocates smaller values of the 
mixing-length parameter over at least the late-A spectral 
range. 

All the models presented here were calculated with no 
overshooting. The baseline grids use l/H = 1.25, but sup- 
plementary grids at l/H = 0.50 were also computed. 
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Figure 1. Grid for solar-abundance models (the lowest gravities 
may be 0.1-0.2 dex larger at higher metallicities). 



2.3 Parameter space 

The grid sampli ng adopted here matches th e basic T E g / log g 
sampling of the ICastelli fc Kuruczl (|2004l ) models; that is, 
from 3.5-13kK at 0.25kK intervals, and 12.0-50.0kK at lkK 
intervals, from logg = 5.0 dex (cgs) to lower values at 0.5- 
dex intervals. This grid density is sufficient for interpolation 
of output products (e.g., linearly in log T e ff/ log g) to be sat- 
isfactory for most applications. Some effort has been put 
into generating models to as low a gravity as possible at a 
given temperature (the effectiv e Eddington limit) , and so 
the present grids extend beyond ICastelli fc Kuruczr s in this 
regard. The sampling is illustrated in Fig. Q] 

Previous work suggests that the Atlas9 turbulent- 
velocity parameter v t has little consequence for limb- 
darkening in the optical regime, and its value (as well as 
its physical significance) is difficult to establish observa- 
tionally in any individual case. Most of the model inten- 
sities presented here have been computed with a canoni- 
cal vt — 2 kms -1 . Solar-abundance grids at vt = and 
4 km s _1 confirm the insensitivity of limb darkening to v t . 

Results of grids of 6571 ne w models are pre- 
sented here. ICastelli fc Kuruczl (|2004l ) h ave calculated mod- 
els us ing solar abundances listed by I Grevesse fc Sauvall 
d 19981). which differ only s lightl y from those given by 
lAsplund. Grevesse fc Sauva 3 (|2005l ). Their grids extend over 
a greater range in scaled solar abundances than those newly 
calculated here, and include additional modifications to 
abundances of alpha-process elements. For completeness, 
intensities have also been calculated from the 8568 atmo- 
spheric structures they provide, along with limb- darkening 
coefficients and other products to match the new grids. 
Table Q] summarizes the grids, where 'Atlas9.C 04' and 
' Atlas9.A10' identify the lCastelli fc Kuruczl (|2004l ) and cur- 
rent grids, respectively. 



3 DATA PRODUCTS 

The principal motivation for the work reported here was to 
provide up-to-date limb-darkening coefficients (particularly 
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Table 1. Summary of ATLAS9 models. The sources for the atmospheric structures are ICastelli fc Kurucj l|2004l ; ATLAS9.C04) and the 
present paper (ATLAS9.A10). In the 'Grid' columns, with entries of the form 'xllvjj', 'xIP indicates the scaling applied to metals from 
the default solar abundance (such that 'mlO' indicates a metallicity [M/H]= —1.0, and l p05' +0.5), while 'vJJ' indicates the turbulent 
velocity used in the model ('JJ' in km s — 1 ). Models at LMC and SMC abundances arc so labelled. The quantity l/H is ratio of mixing 
length to pressure scale height. The number of models in a given grid, jV, is greater for the ATLAS9.A10 set because these generally 
extend to somewhat lower gravities than the ATLAS9.C04 calculations, at all 76 grid temperatures. 



Structure source 


Grid 


l/H 


N 


Grid 


l/H 


N 


ATLAS9.C04 


m25v02 


1.25 


476 


m40av02 


1.25 


476 


ATLAS9.C04 


m20v02 


1.25 


476 


m25av02 


1.25 


476 


ATLAS9.C04 


ml5v02 


1.25 


476 


m20av02 


1.25 


476 


Atlas9.C04 


ml0v02 


1.25 


476 


ml5av02 


1.25 


476 


ATLAS9.C04 


m05v02 


1.25 


476 


ml0av02 


1.25 


476 


ATLAS9.C04 


pOOvOO 


1.25 


476 


m05av02 


1.25 


476 


ATLAS9.C04 


p00v02 


1.25 


476 


p00av02 


1.25 


476 


Atlas9.C04 


p02v02 


1.25 


476 


p02av02 


1.25 


476 


Atlas9.C04 


p05v02 


1.25 


476 


p05av02 


1.25 


476 


Atlas9.A10 


POOvOO 


1.25 


554 


LMCv02 


1.25 


554 


Atlas9.A10 


P00v02 


1.25 


554 


LMCv02r 


0.50 


552 


Atlas9.A10 


P00v02r 


0.50 


552 


SMCv02 


1.25 


554 


Atlas9.A10 


P00v04 


1.25 


548 


SMCv02r 


0.50 


552 


Atlas9.A10 


P03v02 


1.25 


543 








Atlas9.A10 


P03v02r 


0.50 


544 








Atlas9.A10 


P05v02 


1.25 


531 








Atlas9.A10 


P05v02r 


0.50 


533 









at MC abundances). Thus the main data products are spe- 
cific intensities (monochromatic radiances) I\((J,), as a func- 
tion of the cosine of the angle between the line of sight 
and the surface normal. These are provided in condensed 
form as broad-band limb-darkening coefficients. 
Physical fluxes, 



been used: 



Fx = 2tt / / A (/-)/- d/_ 
Jo 



(1) 



(where H\ is the Eddington flux), are required as a check on 
the accuracy of intensity integrations, and in any case repre- 
sent only a minor additional computational overhead. Model 
fluxes are therefore also provided, both at the 1221 wave- 
lengths used in the structure calculations (9lA-160/im), and 
in broad-band ('synthetic photometry') form. Other data 
products (such as model structures and higher-resolution 
intensities) will be made freely available on request. 



4 LIMB DARKENING 
4.1 Passbands 

Modern detectors working in the optical are essentially 
photon-counting in nature. For such detectors broad-band 
model intensities (and, mutatis mutandis, fluxes) may be 
calculated as 



hiii) 



f x i x {vL)<l> x (i)\d\ 

/ A 0A(i)AdA ' 



(2) 



(e.g.. lBesselill2005h where (f> x (i) is the response function of 
passband i in a given photometric system. Historically, the 
equivalent formalism for energy-integrating detectors has 



Ufa) 



fx h(n)(j>x(i)<i\ 
IxMi)dX- 



(■■>>) 



Results are provided for photon-counting detectors for 
all photometric systems considered here, and addition- 
ally for energy-integrating detectors for the Johnson and 
Stromgren passbands. The adopted sources for response 
functions are: 

• UBVRI (Johnson-Cousins system) : IrBesselll (Il990l) 



1 198c 



JHKLL'M (Johnson-Glass system): iBessell fc Brett] 



uvb y (Stromgren s ystem ) : ICrawford fc Barnes! (|l970l : 



filters), IStubbs et al 



l systerr 
] |2007l : 



at mosph e ric t ransmission) , 



iHamamatsul (|l999l ; S 4 response), I Alien! (Il976l, p. 108; alu- 
minium reflectivity), IPraezisions Glas fc Qptikl (BK7 glass 
transmission) 

• ZYJHK (WFCAM system) : iHewett et al.l (l200f]|) 

• H P! Bt, Vt (Hipparcos/Tycho system): iBesselll (|2000h 

• ugriz (Sloan system): on-line, 

http : //www . sdss3 . org/ instruments/camera. php 

• Kepler: on-line, 

http : //keplergo . arc .nasa.gov/kepler_response_hiresl . txt 

Simple geometric treatements of mutual irradiation (the 
'reflection effect') in binary-star systems r equire the an gu- 
lar dependence of the bolometric intensity (|Wilsonlll990] ) , so 
bolometric limb-darkening coefficients have also been calcu- 
lated (<t>x = 1). 

2 The Stromgren system is close to being filter defined, but at- 
mospheric transmission, 1P21-S4 photomultiplicr sensitivity, and 
BK7 glass transmission all make roughly equal inroads into the 
short- wavelength edge of the u passband. 
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4.2 Analytical characterization 

Specific intensities are traditionally, and conveniently, char- 
acterized by simple functional forms, of which the most ven- 
erable is the linear limb-darkening law 



= i - «(i - aO, 



(4) 



the analytical s olution for a source function that is linear in 
optical depth r jSchwarzschildll 19061 iMilneil 192 it u = 0.6 for 
a grey atmosphere). 

More-realistic models do not have analytical solutions, 
and additional terms are required for a faithful func tional 
representation of actual limb darkening. iKopaJl (|l949T ) took 
a series-expansion approach, and found the differences be- 
tween quadratic and quartic approximations to be neg- 
ligible at the level of accuracy with which he was con- 
cerned. Two-coefficient parameterizations subsequently be- 
came standard; in particular, a quadratic law of the form 

1(^/1(1) = l-a(l- M )-&(l- M ) 2 (5) 

was widely adopted in the modern computational era (e 
iManduca. Bell fc Gustafssonlll977l : IWade fc Rucinskilll98. 
although a cubic version, 

/( M )//(l) = l-a(l- M )-6(i-M) 3 , 

was advocated by Ivan 't Veer! l|l960l ). An alternative loga- 
rithmic law, 

7( A1 )/J(l) = l-a(l-p)-6(/xln/x) (6) 

was proposed by iKlinglesmith fc Sobieskil (|l97Cf ). and a 
square-root law, 



■ a(l - fj) - 6(1 - y/JI), 



bv lDi'az-Cordoves fc Gimenej (|l992f ). 

For modern work, these limb-darkening laws suffer from 
various limitations. The obvious way to develop a single an- 
alytical expression that more accurately reproduces numer- 
ical limb-darkening results across a wide parameter s pace is 
simply to use more terms, a notion implemented by IClaretl 
who introduced a four-coefficient fit in powers of 



4t 

J(/0/I(l) = l-£a»(l 



/'■' 



i/2 



(8) 



Figure [5] illustrates the differing degrees to which the 
foregoing common functional forms reproduce the intensities 
computed here. While quantitative results depend on the 
chosen normalization and adopted figure of merit, the qual- 
itative conclusions are robust: a linear limb-darkening law 
is a rather poor representation of detailed numerical results; 
a quadratic law fares little better; but the four-coefficient 
law affords an improvement of up to an order of magnitude 
in precision over any standard two-coefficient option. While 
there has been some discussion in the literature as to which 
is the most appropriate limb-darkening 'law' to adopt for 
light-curve analyses, and in particular under circumstances 
which permit e mpirical determin ation of limb-darkening co- 
efficients (e.g.. [Southwortb"ll2008f ). there seems little justifi- 
cation for adopting anything other than a high-order fit such 
as eqtn. (JSj) if the principal purpose is to secure an accurate 
representation of model-atmosphere results. 




4— coefficent 



[M/H] = + 0.0 v t =2 km/s 
log(g)=4.5 I/H=1.25 
R band 
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Effective Temperature (kK) 

Figure 2. Comparison of different parameterizations of limb 
darkening (equations [4j|8j) . The quantity plotted is log 10 of the 
maximum absolute difference between input and parameterized 
values of I(fi)/I(l). 
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Figure 3. Differences, flux minus integrated intensity (expressed 
as a percentage of flux) for several limb-darkening laws (see Sec- 
tionl4~3ll. 



4.3 Fitting procedure 

For any given analytical form, there remains a choice of nu- 
merical technique for determining the limb-darkening coeffi- 
cients. The most direct approach is simple least-squares de- 
termination of coefficients from functional fits to computed 
values (with 1(1), the fitted value of 1(1), as an op- 
tional additional free parameter). A r elated approach, the 
r-integration method dHevrovskvl l 2007h . appears to offer no 
significant benefits (|Claretll2008l ). 

A criticism that can be levelled against the least-squares 
method is that subsequent integration of parameterized in- 
tensities may not accurately recover the flux, as it should 
(eqtn. [!}. Alternative formulations in w hich which flux is 
explicitly conserved can be devised (e.g. JWade fc Rucinskil 
Il985l ; lvan Hammeill99 j ). giving rise to so-called flux conser- 
vation methods (cf., e.g., IClaretl 120081 . for a recent discus- 
sion )^_SjJ£h_jnethodsreqtm further ad hoc constraint; 
e.g.. IWade fc Rucinskil (|l985l ) set 1(1) = 1(1) for the linear 
model. 

For the most part, straightforward least-squares fitting 
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has been used here, to intensities evaluated at 20 angles[f| 
Subsequent analytical integration of the four-coefficient in- 
tensity model then reproduces directly computed fluxes to 
better than 0.1%. When fit in this manner, all two-coefficient 
models conserve flux to within ~1%, while the linear model 
again fares poorly (Fig. 0. 



5 PHOTOMETRIC MATTERS 

Broad-band fluxes have potential utility in a variety of ap- 
plications (effective-temperature determinations, population 
synthesis, etc.). Thus while it is not the intention here to 
conduct a detailed confrontation of the new model fluxes 
with observations, a basic transformation of the synthetic 
photometry to observed magnitude systems may be of inter- 
est. Modern systems generally have dedicated programmes 
of flux calibration, so the focus here is on the Johnson and 
Stromgren systems. 



5.1 Johnson system 

5.1.1 Flux zero-point 

Full normalization of the synthetic photometry can be di- 
vided into two parts: absolute flux calibration, and differ- 
ential (colour) calibration. The former aspect is tradition- 
ally bound to the absolute flu x calibration of V ega, which 
has be en addressed in detail bv lHavesl (| 19851 ) and lMegessierl 
|l995l ). using primary measurements mostly made in the 
1970s. Their key results can be summarized as pseudo- 
monochromatic 5556-A fluxes: 



/(5556) = 3.44 x 10~ 
3.46 x 10~ 



erg cm 2 s 1 J 
W m~ 2 nm _1 



respectively. These results can be used to provide a zero- 
point calibration for the present grids by using a model 
matching Vega, and its observed V-band magnitude. 

Both 'the' V-band magnitude, and the choice of model, 
are s ubject to uncertainty ; widely adopted values are V = 
0.03 (| Johnson et all 1 19661 ), and r ? ff = 9550K, logg = 3.95 , 
[M/H] = -0.5, v t = 2 kms -1 (|Castelli fe Kurud fl994[ ) . 
Subsequent a nalyses have supported this effective tem- 
perat ure (e.g. ICiardi et al.l 1200 ll : iHill. Gulliver fc Adelmanl 
2010), but the HST calibration programme uses a model 
with T ef f = 9400 K, lo g. q = 3.90, [M/H] = -0.5, v t = 
kms" 1 (|Bohlin|[2007l ). Workers are therefore faced with 
choices in observed flux, observed magnitude, and adopted 
model, as well as in whether to adopt energy-integrating or 
photon-counting magnitudes. The situation is further com- 
plicated by the fac t that Vega is a rapid rotator viewed 
pole-on dGravlll98g|; [Gulliver et - al.il 199 if : lAufdenberg et al] 
120061 : IHill. Gulliver fc Adelmanl |2010T ). so that in principle 
no single-temperature model can give an accurate represen- 
tation of the entire flux distribution (although for magnitude 
zero-point calibration this is of minor importance, because 
the models are only used to scale between 5556-A and V- 
band fluxes). 



3 For comparison purposes, the linear law has been modelled with 
and without 1(1) as an optimized parameter, and with explicit 
flux conservation; Appendix FBI 



The baseline adopted here i s the Hayes calibra- 
tion, 9.55kK model, Ijohnson et al.l (1 19661 ) V, and photon- 
counting magnitudes. For this system the V-band flux f(V) 
is 1.031 x /(5556), yielding 



V = — 2.51og 10 f(V) — 21.096 

+ C1+C2 + c 3 + c 4 



(9) 



where f(V) is in erg cm 2 s 1 A 1 and the optional correc- 
tion factors are 

Ci = -0.0004 for the T ef f = 9400 K m odel; 

C 2 = +0.0063 for the lMegessierl (| 19951 ) /(5556) flux; 

more generally, C2 = 2.51og 10 [/(5556)//(Hayes)]; 
Cs — +0.0135 for energy-integrated photometry; and 



C4 = Vvc: 



0.03. 



5.1.2 Colour-index zero-points 

Differential normalization relies on relating the observed 
colours of a standard star, or stars, to model fluxes. A model 
colour index is formed from fluxes in two passbands PI and 
P2by 



(P1-P2)* 



-2.51og 10 
-2.51og 10 



/GP2) 



+ {Cpi — Cp-i) 

+ Z P1 2 (10) 



where /(P1),/(P2) are model broad-band fluxes and 
Cpi, Cp2, Zpvi are constants. Equating the right-hand side 
of eqtn. (|10[) with the observed (reddening-frce) colour index 
yields the required normalizing constant Zp\2. 

Vega has traditionally been adopted as the primary cali- 
brating star, with increasinging u se of Sirius as a supp lemen- 
tary standard since the work of I Cohen et all (|l992h . Once 
again, however, observed photometry is subject to both sta- 
tistical uncertainties and deterministic changes (e.g., one 
might choose to define Vega to have colour indexes of zero, 
or to average values from an ensemble of standard stars), 
and the appropriate model parameters are also moo t . Here 
the photometry compiled bv lBessell. Castelli fc Plezl (|l998l ) 
is adopted, as are the Sirius model parameters due to Ku- 
rucz which they quote: T cff = 9.85k K, logff = 4.25 , [M/H] 
= +0.5, Vt = 2 km s™ (see also ICastellil 1 19991 ). These 
parameters are in good agreement with other recent de- 
terminations Je.g.. lHill fc Landstreetlll993l ; Ivan Noortlll998l ; 
iQiu et alJlioOll ). 

The resulting normalizing constants are listed in Ta- 
ble [2] Agreement between the Vega and Sirius calibrations 
for the normalizing constants inferred from the specific set 
of observed colours adopted here is better for the 9.55kK 
Vega model, and the adopted Zp\2 values are the average of 
this model and the Sirius results. 



5.1.3 Bolometric corrections 

From the basic definitions of bolometric correction and mag- 
nitude, it is straightforward to show that the bolometric 
correction in some passband P for a model at effective tem- 
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Table 2. Johnson colour- index normalizations. Column 2 gives the colour zero-point computed from the Vega 9.55kK model, with 
photometry adopted from [Bessell. Castelli fc Pled jl998l : e.g., observed (V — R) = —0.009). Column 3 gives the equivalent results for 
Sirius. The recommended zero-points, Zpi2, are the average of the previous two columns (with corrections for alternative photometries 
dropped); synthetic colours may then be computed from model fluxes by using cqtn. (1 10K . 

Ce is the correction to be added to columns 2-4 for energy-integrating detectors, and Cg.4 is the c orrection to be added to column 2 to 
obtain results for the 9.40kK Vega model. The Ux and Bx passbands are those recommended bv lBesselll l ll990l) for use in computing- 
synthetic colours to be compared with observed (U — B) colours. 



Colour 




Vega 9.40kK 




Sirius 9.85kK 




Zpi2 


C E 


C9.4 


(Ux 
(B 


-B x ) 
-V) 


-0.422 
+0.609 


+ [(U-B)- 
+ [(B -V)- 


0.000] 
0.000] 


-0.418 
+0.608 


+ [(u 

+ [(B 


-B) + 
-V) + 


0.045] 
0.010] 


-0.420 
+0.608 


+0.031 
+0.000 


+0.020 
+0.007 


(V 
(V 
(V 


-R) 

-I) 
-K) 


+0.568 
+1.260 
+4.924 


+ [(V - R) + 

+ [<y-i) + 

+ [(V — K) - 


0.009] 
0.005] 
0.020] 


+0.575 
+1.276 
+4.918 


+ VY 

+ [(V 
+ [(V 


-R) + 

-I) + 
-iv)4 


0.010] 
0.016] 
0.061] 


+0.572 
+1.268 
+4.921 


+0.018 
-0.002 
-0.002 


+0.005 
+0.010 
+0.026 


(J 
(H 
(K 


-K) 
-K) 
-L) 


+2.261 
+ 1.156 
+1.872 


+ [(J-K)- 
+ [(if - K)} 
+ [{K — L) — 


0.010] 
0.010] 


+2.249 
+1.152 
+1.872 


+ [(J 

+ [{H 
+ [{K 


-K) + 
-K) 4 
— L) — 


0.018] 
- 0.009] 
0.003] 


+2.255 
+1.152 
+1.872 


-0.002 
+0.002 
-0.001 


+0.006 
+0.002 
+0.002 



perature T B g is given by 

BC(P) =M Bo i,0- 2.5 log 10 



I K (3 



2.5 log 10 T e 4 ff + 2.5 log 10 (f P ) + C P 



(11) 



where dio is 10 pc (expressed in the same units as the solar 
radius, Rq) and Cp is the magnitude- normalizing constant, 
which can be found for the Johnson passbands from data 
summarized in the preceding two subsections|j] 

Bolometric corrections in the main Johnson bands are 
included for each model in the on-line listings, adopting 
MboI,© = +4.74. As a check on results, a solar-like model 
was calculated with T cH — 5778K, logp = 4.44 (cgs), 
v t = 2 km s" 1 , i/H = 1.25, [M/H] = 0.0; for this model, the 
(photon-counting) l^-band bolometric correction is found to 
be BC(V) = —0.075, in ex cellent agreement with the em - 
pirical solar value of -0.07 (|Bessell. Castelli fc Plezlll99Sl) . 



5.2 A Stromgren calibration 

Flux calibration is readily carried out for the Stromgren uvby 
system in the same way as for the Johnson system. A po- 
tential complication is that for the intermediate-width pass- 
bands of this system, individual spectral features may be sig- 
nificant; in particular, H5 (410nm) falls in the middle of the 
v band. To check that the standard 20-A spectral sampling 
is adequate in these circumstances, high-resolution synthetic 
spectra were computed from the Vega models. Stromgren- 
band fluxes from these synthetic spectra agree with the stan- 
dard results to better than 0.1% for u, b, and y, and are only 
0.5% smaller at v. Since the Balmer lines reach their great- 
est strength at early-A spectral types, the standard sampling 
appears to be satisfactory. 

Results are summarized in Table [3] Stromgren colours 



4 E.g., C B = (C B - C v ) + C V = +0.608 - 
where numerical values for photon-counting 
Tableland eqtn. ©. 



21.096 = -20.488, 
photometry are from 



can be generated from the models with 

7W 



(b-y) 



m 1 



ci = 



-2.51og 10 
-2.5 log 10 
-2.51og 10 



+ Z, 



m 

f(v) x f(y) 



by 



P(b) 
/(«) x f(b) 



+ Zm\ 
+ Z cl 



where f(p) is the model flux in passband p and the zero- 
poi nts Z axe the last t hree e ntries in column 2 of Table [3] 
iFabregat fc Reid 

CSSi) and iGravl { 19981 ) have previ- 
ously provided Stromgren flux calibrations which are the 
results of convolving th e uvby passbands with o bserve d spec- 
tmphotometryQ from iGlusfmeva et al.l (ll992Tl and iTavlorl 
l|l984h . respectively (both nominally on the [Hayes 1985| sys- 
tem). Comparison with their results, expressed as a zero- 
magnitude flux in units of 10 erg cm" - 2 s- 1 A' 1 for consis- 
tency with their presentations, shows generally good agree- 
ment: 



Band 


FR 


Gray 


Atlas9.A10 


u 


11.80 


11.72 


12.20 


V 


8.69 


(8.66) 


8.73 


b 


5.84 


5.89 


5.90 


11 


3.69 


3.73 


3.72 



The ~4% difference at u between spectrophotometric and 
model results may be a consequence of uncertainties in 
observed Vega magnitudes, which range over 10% in this 
passband (|Perrvll969l;lBarrvll969l ; |Johansen fc Gvldenkernd 



ll97d ; ICrawford fc Barnes!ll970l ); e.g., using Sirius as the cal- 
ibrator yields a zero-magnitude u-band flux of 11.79 x 10~ 9 



erg cm 2 s 1 A 



6 DISCUSSION 



Much attention has been paid in the literature as to which 
is the best analytical representation of limb darkening, and 

5 The Gray i;-band calibration is model-based. 
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Table 3. Summary of Stromgren photometric calibrations (see Sect ion |5. 211, The second colum n gives flux-calibratio n and c olour zero- 
points for photon-counting detectors, using a T c g = 9.55kK model. ICrawford &: Barnes! jl970f ) photometry and the iHavesl j 19851 ) flux 
calibration for Vega. 

C\ is the correction to be added to the values in column 2 if adopting other sources of photometry; 

Ci is the correction to be added for ene rgy-integrating d etectors; 

Cz is the correction to be added for the lMegessierl l|l995h flux calibration. 
Column 6 gives the corrections to be added for a T e ff = 9.40kK Vega model, while the final column gives the corrections for Sirius 
photometry (baseline data from lCrawford. Barnes fc G olson 197(J) and a 9.85kK model, where AV = Vye ga — 0.03. 



Band/ 


Vega 


Ci 


C 2 


C-i 


Vega 




Sirius 


Colour 


9.55kK 








9.40kK 




9.85kK 


U 


-19.784 


+(u - 1.445) 


+0.001 


+0.0063 


-0.045 


+AV - 


0.037+ (u - 1.298) 


V 


-20.147 


+(v - 0.195) 


+0.000 


+0.0063 


-0.016 


+AV - 


0.019 + (v- 0.166) 


b 


-20.573 


+(6 - 0.034) 


+0.001 


+0.0063 


-0.011 


+AV + 


0.012 + (6- 0.036) 




-21.074 


+(V - 0.03) 


+0.002 


+0.0063 


-0.006 


+AV - 


0.002 


(b-y) 


-0.500 


+0.004 -(b-y) 


+0.000 




+0.004 


-0.015 


+ [0.006 - (b - y)\ 


m,\ 


+0.074 


+0.157 - mi 


+0.001 




+0.001 


+0.047 


+ [0.124 - mi] 


CI 


+0.064 


+1.089 - ci 


-0.002 




+0.023 


-0.015 


+ [1.002 - ci] 



[M/H] = + 0.0 v.=2 km/s 
R band l/E-1.25 - 




5.0 6.0 7.0 8.0 9.0 

Effective Temperature (kK) 

Figure 4. Maximimum differences in I (fi)/I(l ) betw een the re- 
sults of the present models and fits by [claret ( 2000]) to earlier 
Kurucz models. 



which is the appropriate numerical technique for evaluating 
limb-darkening coefficients (see discussion and references in 
Sections l4.2ll4.3[) . While these are important details, the use 
of a high-order parameterization, which provides satisfac- 
tory results under most circumstances, renders such details 
academic. 

However, the larger importance of the underlying model 
atmospheres has not generally been emphasized, perhaps 
because most analysts have used 'off-the-shelf model- 
atmos phere results. The new o pacity sources introduced 
in the ICastelli fc Kuruczl (|2004l ) models show up directly 
in emergent fluxes and intensities, and the changes in 
treatment of convection also change the structures. These 
changes are particularly important at T c g < 8kK; for mod- 
els at T e ff > lOkK the differences bet ween the curre nt mod- 
els, th ose of lCastelli fc Kuruczl (|2004l ), and the older [Kuruczl 
grids, are small, and result largely from revisions to 
adopted abundances. 

As far as specific intensities are concerned, differences 
between the newer calculations and older results are illus- 
trated by a comparison between the extensive listings pro- 
vided by IClaretJ (|2000l ) and the Atlas9.A10 models pre- 
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0.8 
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— • — &/H 1.25 - 
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_ v t =2 km/s 




1 1 — 1 — — 1 — 1 1 — — + 

_ 3.5kK, 4.5 
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Claret (2000) - 





0.0 



0.2 



0.4 0.6 

ix = cos(7) 



0.8 



1.0 



Figure 5. Normalized i?-band intensities from the present calcu- 
lations for two T e ff/logg values. Filled circles are calculated in- 
tensities, andcontinuous lines the corresponding four-coefficient 
fits. The lClareq l|2000h results characterize earlier Atlas9 models. 



sented here (Fig. UJ. The changes can be much larger than 
any that arise through different fitting techniques, provided 
that the adopted functional form is adequate (Fig. [SJ. 

We might expect further developments in Atlas mod- 
els in future, but there is reason to hope that consequences 
may be modest in scale, at least in the short term; using 
the latest version of his opacity-sampling Atlas12 code, 
Kurucz (personal communication) has kindly computed a 
solar-abundance model at T e ff = 7.5kK, logg = 4.0, ft = 
2 km s . Differences between old and new Atlas9 models 
are relatively large for these parameters, but further changes 
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in the Atlas12 products are negligible for the present pur- 
poses. 

Finally, the limitations of the present work should be 
recognised. First, although the Kurucz/Castelli ODFs used 
here represent a significant improvement over older compi- 
lations, they are still inevitably incomplete. In particular, 
there are a number of indications that Atlas9 models be- 
gin to show significant disagreeme nt with observations for 
effect i ve temperatures <4 kK (e.g., iBessell. Castelli fe Plezl 
1 19981 ; iBert one et ail 120081 ). presumably because of the in- 
creasing importance of missing molecular opacities. Models 
at the coolest effective temperatures should be used with 
caution. 

Secondly, the synthetic photometry can be expected 
to show not only the successes i n reproducing observa- 
tions repor ted for similar m odels bv lBessell. Castelli fe Plea 
Jl99ST) and ICastellil ( 1999), but also the failures; note that 
IBessell, Castelli fc Pled (|l99Sl . cf. their Appendix E.3.1) find 
that their modelled (U — B) colours (i.e., using Ux — Bx 
passbands) required an ad hoc scaling of 0.96 x to bring 
them into agreement with observations. 

Thirdly, it must be owned that the new models pre- 
sented here fare no better than old ones in matching 
the handful of apparently reliabl e empirical determina- 
tions of ste llar limb darkening (e.g. . lHevrovskvll2007l ; Iciaretl 
l2008l , 120091 ). although such determinations are usually lim- 
ited to estimates of coefficients for linear or, occasion- 
ally, quadratic limb-darkening laws for stars other than 
the Sun (and such laws are poor representations of model 
results). Common simplifying model-atmosphere approxi- 
mations for scattering phase functions will be a factor in 
this; although such approximations may be appropriate for 
structure and integrated-flux calculations, they have obvi- 
ous weaknesses for evaluating intensities (or polarizations). 
Addressing these issues is a matter for future study. 



7 RECOMMENDATIONS 

Practically all modern studies in the optical region use 
photon-counting detectors, and so models appropriate to 
that context will normally be used. The investigator is still 
faced with choices in model convection, microturbulence, 
and abundance (Table [I}, and in limb-darkening parameter- 
ization (equations rjHH])- If, as is frequently the case, there 
are no results from detailed analyses to suggest otherwise, 
it's reasonable to adopt the appropriate global abundance 
system (solar, LMC, or SMC), v t = 2 km s~\ t/H = 1.25. 
The careful investigator might well conduct sensitivity tests 
to examine the consequences of other parameters (including 
varying T E g and log<?); in particular, for effective tempera- 
tures less than ~8kK, it would be prudent to examine the 
effects of models computed at l/H — 0.50. 

Other than for comparison with empirically-determined 
limb-darkening laws (normally linear, but increasingly for 
quadratic coefficients), there is little merit in adopting any- 
thing other than a four-coefficient limb-darkening parame- 
terization (equation [8| if seeking an accurate representation 
of model-atmosphere intensities. The computational cost 
compared to simpler formulae is trivial in anything other 
than Monte-Carlo calculations, and even there the overhead 
should be entirely tolerable. 
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Table Al. Adopted LMC and SMC abundances, by number, on a scale where the hydrogen abundance is 12.0 dex. 



Element LMC SMC Element LMC SMC Element LMC SMC Element LMC SMC 
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He 


10.88 


10.85 


3 


Li 


0.75 


0.45 
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JJe 


1.08 


0.78 
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jj 


2.40 


2.10 
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7.73 


7.37 
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6.90 


6.50 
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8.35 


7.98 


Q 


w 
I 




o.yo 


10 


Ne 


7.60 


7.20 


11 


Na 


6.97 


5.78 


12 


Mg 


7.06 


6.72 


13 


Al 


6.07 


6.27 


14 


Si 


7.19 


6.79 


15 


I' 


5.06 


4.76 


16 
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6.70 


6.30 


17 


CI 


5.01 


4.95 


18 


Ar 


6.20 


5.90 


19 
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4.78 


4.48 


20 


Ca 


5.83 


5.63 


21 


Sc 


2.62 


2.31 


22 


Ti 


4.63 


4.31 


23 
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4.10 


3.56 



21 


Cr 


5.43 


5.06 


25 


Mn 


5.06 


4.88 


26 


Fc 


7.23 


6.93 


97 
z < 


Prt 


4 62 


4 32 


28 


Ni 


6.02 


5.82 


29 


Cu 


3.91 


3.63 


30 


Zn 


4.20 


4.00 




V.T<t 
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APPENDIX A: ADOPTED MAGELLANIC-CLOUD ABUNDANCES 

Baseline solar abundances were taken from Asplund. Grevesse fc Sauvall (20051) . Cloud abunda nces were adopted , in order 
of pre ferenc e, from: FLAMES result s (|Hunter et al.ll2007l . iTrundle et alj|2007l , as summarized bv lEvans et al.ll2008T ); ICarnettJ 
(1999j); and iRussell fc Dopital (1 19921 . applying their differential values to Asplund solar abundances). In each case, samples 
of early-type stars or H II regions were analysed, so the abundances are representative of Popul ation I. For elements no t 
studied in these sources, solar abundances were used, adjusted by the median metal offsets from IRussell fc Dopital (|l992f ): 
—0.3 (LMC) and —0.6 dex (SMC). These global adjustments are in adequate agreement w ith more-recent discussions (e.g., 
iHill. Andrievskv fc Spitell 19951 ; [Hu1|199T| ; iLuck et al.lll998l ; IVennlll999l ; iMokiem et aT1l2007f l. 
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APPENDIX B: ON-LINE MATERIAL 

New opacity-distribution functions using the lAsplund. Grevesse fc Sauvall (|2005l ) solar abundances are provided on-line, at 
metallicities [M/H] = -0.5, +0.0, +0.2, +0.3, +0.4, +0.5, together with the LMC and SMC ODFs. In each case, ODFs are given 
sampled at the standard 'big' and 'little' wavelength sampling intervals required for use with Atlas9, and for microturbulent 
velocities of 0, 1, 2, 4, and 8 kms -1 , together with the corresponding Rosseland-opacity files. These are binary files in the 
format required by Atlas9. 

Limb-darkening and flux d ata products are organiz ed into directories corresponding to the grids listed in Table [T] (for 
both new models, and those of ICastelli fc Kurucj 12004 ) . Within a given directory, there are three plain-text files for each 
model, with names of the form t05000g40 . <ext> where the base identifies T e g and logg, and the extension identifies the file 
content. 

Files with .fix extensions contain listings of physical fluxes at 1221 wavelengths, from 90A to 160^im (with 20-A sampling 
through the optical). Fluxes are tabulated in erg cm~ 2 s _1 A -1 as a function of wavelength (in A). 

Limb-darkening coefficients are provided in files with extensions .ucP (for photon-counting detectors; Section [4.1[) and 
.ucE (for energy-integrating detectors). The .ucE files list broad-band limb-darkening coefficients for Stromgren and Johnson 
systems (including Ux and Bx passbands for computing synthetic (U — B) colours), while the .ucP files give broad-band 
limb-darkening coefficients for all 30 passbands listed in Section T4. II plus [energy- integrated] bolometric coefficients. The data 
given for each passband may be illustrated by V-band results for a Vega model (T e ff = 9.55kK, logg = 3.95, [M/H]= —0.5, 
v t = 2 kms" 1 , £/H = 1.25): 
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Here the first line lists 

(i) the passband, 

(ii) effective wavelength (in A), 

(iii) the broad-band physical flux /(A) (in erg cm -2 s _1 A -1 ), 

(iv) the broad-band surface-normal intensity (in erg cm -2 s _1 A -1 sr _1 ), and 

(v) the bolometric correction (for Johnson passbands only). 

Subsequent lines list coefficients for the limb-darkening laws summarized in Section 14.21 the final two numbers in each row 
are the r.m.s. and maximum differences between input and modelled values of 

The 'linear - FC and 'quad - FC results use flux conservation in place of least squar es, with the fitted valu e 1(1) 
anchored to the model 1(1) in the first case, and additionally 1(0.1) to 1(0.1) in the second fcf. lWade fc Rucinski|[l985l ) . The 
'linear - 1' listing is for the standard (single-parameter) least-squares fit to eqtn. Q, while the coefficients in the 'linear 
- 2' listing are for a model with 7(1) free, i.e., 

I( l i)/I(l) = a-b(l- l x). 
This paper has been typeset from a TpX/ F/IpX file prepared by the author. 



